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Abstract 

The  focus  of  this  research  is  an  investigation  of  energy  band  anisotropy  in 
simple  quantum  well  structures.  This  anisotropy  results  from  the  asymmetry  of 
the  periodic  potential  within  the  crystal  lattice.  For  sufficiently  high  doping  levels, 
band  structure  anisotropy  is  expected  to  play  an  important  role  in  the  evaluation  of 
the  electronic  and  optical  properties  of  the  quantum  well  structures.  The  analysis 
uses  a  model  based  on  a  6x6  Luttinger-Kohn  kp  approach  for  bulk  material  valence 
band  structure  together  with  the  Envelope  Function  approximation.  The  model  is 
used  to  analyze  Si/Sii-xGe*,  GaAs/AlxGai-*  As,  and  GaAs/InxGai_xAs  quantum 
wells  for  the  [110]  and  the  [001]  growth  directions.  The  resulting  bandstructures 
show  significant  anisotropy  for  materials  grown  in  both  the  [110]  and  [001]  growth 
directions.  In  all  cases  the  materials  grown  in  the  [110]  direction  show  a  more 
pronounced  anisotropy  than  the  materials  grown  in  the  [001]  directions.  For  the 
[001]  growth  directions,  the  bandstructures  were  effectively  isotropic  for  values  of  ky 
less  than  0.4  Angstrom -1  for  Si/Si7oGe3o,  0.6  Angstrom _1  for  GaAs/Al3oGa7oAs, 
and  0.5  Angstrom _1  for  GaAs/In3oGa7oAs. 
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Band  Structure  Anisotropy  in  Semiconductor  Quantum  Wells 


I.  Introduction 

Since  the  inception  of  the  semiconductor  quantum  well  in  1970  by  Esaki  and 
Tsu  (12)  and  the  later  advances  in  Molecular  Beam  Epitaxy  (MBE)  and  Metal- 
Organic  Chemical  Vapor  Deposition  (MOCVD),  quantum  well  devices  have  proven 
to  be  very  useful  in  the  field  of  optoelectronics.  Two  significant  devices  in  which 
these  structures  are  used  are  the  laser  diode  and  the  infrared  photodetector  (15). 

Of  particular  interest  to  the  Air  Force  is  the  use  of  these  bandgap-engineered 
materials  in  infrared  photodetectors  in  the  wavelength  range  corresponding  to  the 
atmospheric  windows  of  3-5 um  and  8-14 /xm  for  use  in  applications  such  as  thermal 
imaging,  guidance,  and  reconnaissance  (2), (16).  Quantum  well  photodetectors  are 
also  well  suited  for  use  in  the  1.3 /xm  and  1.5//m  range  which  makes  them  useful  for 
fiber-optic  communications  (15). 

Quantum  devices  offer  many  advantages  over  normal  photodetectors.  First 
they  have  a  great  flexibility  in  design  which  allows  tailoring  to  a  specific  operating 
wavelength.  As  explained  in  the  next  section,  parameters  such  as  well  width  and 
alloy  composition  can  be  modified  to  change  the  quantum  well  energy  levels  which 
in  turn  change  the  optical  response.  Second,  for  inter-subband  transitions  quantum 
well  devices  have  a  narrow  pass-band  spectral  response.  This  is  a  result  of  the  well 
possessing  a  limited  number  of  energy  levels  between  which  transitions  can  occur 
(15).  Finally,  by  using  materials  such  as  Si/SiGe,  operating  temperatures  can  be 
raised  to  approximately  75K  from  that  required  by  bulk  material  detectors,  which 
is  in  the  region  of  30K  (2). 


1 


Figure  1  Growth  layers  in  an  AlGaAs/GaAs  quantum  well  (19). 

1.1  Quantum  Well  Structures 

Analysis  of  the  optical  properties  for  a  material  typically  requires  knowledge 
of  the  material’s  bandstructure.  The  bandstructure  can  be  modified  by  creating 
heterostructures  which  are  discontinuities  in  the  crystal  created  by  placing  two 
different  materials  in  contact  (10).  The  particular  heterostructure  studied  in  this 
research  is  the  quantum  well.  The  bandstructure  of  the  quantum  well  in  most  cases 
is  very  different  from  that  found  in  the  bulk  material. 

Semiconductor  quantum  wells  are  created  by  alternating  semiconductor  ma¬ 
terials  with  differing  bandgaps.  Through  approaches  like  MBE,  very  thin  layers 
are  created  which  approach  the  deBroglie  wavelength  of  the  electron.  This  allows 
for  quantum  confinement  effects  to  become  important.  The  growth  layers  for  a 
GaAs/AlGaAs  multiple  quantum  well  are  shown  in  figure  1. 

In  crystal  structures  the  periodicity  created  by  the  particular  atomic  arrange¬ 
ment  creates  energy  bands  that  determine  the  electronic  properties  of  the  material. 
This  bandstructure  identifies  the  allowed  energy  eigenstates  for  a  particular  value 
of  the  wave  vector  k.  The  bandstructure  is  commonly  displayed  as  a  function 
of  direction  in  fc-space  from  one  symmetry  point  to  another  within  the  Brillouin 
zone.  The  Brillouin  zone  is  defined  as  a  single  unit  cell  of  the  reciprocal  lattice. 
Examples  of  some  symmetry  points  and  the  first  Brillouin  zone  for  the  zincblende 
crystal  structure  are  shown  in  figure  2. 
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Figure  2  First  Brillouin  zone  and  symmetry  points  for  zincblende  and  diamond 
structure  (22). 

In  a  quantum  well,  the  bandstructure  can  be  modified  to  optimize  the  optical 
properties  of  the  device.  Well  width,  growth  direction,  and  alloy  composition  can 
be  varied  to  maximize  absorption  for  a  particular  wavelength.  Also,  by  creating 
multiple  wells  of  different  widths  and/or  compositions  a  detector  can  respond  to 
multiple  wavelengths  (2). 

Quantum  well  devices  are  usually  comprised  of  many  quantum  wells  and, 
depending  on  the  spacing  between  the  individual  wells,  are  divided  into  two  cate¬ 
gories:  superlattices  or  multiple  quantum  well  structures.  If  the  spacing  between 
the  individual  wells  is  small  enough  such  that  the  wavefunctions  of  the  confined 
carriers  interact  with  the  wavefunctions  of  those  carriers  in  the  neighboring  wells, 
then  a  superlattice  is  formed.  If  the  spacing  between  wells  is  great  enough  such  that 
there  is  essentially  no  interaction  between  wavefunctions  of  neighboring  wells,  then 
the  structure  is  considered  a  multiple  quantum  well.  The  two  types  of  structures 
are  shown  in  figure  3. 

In  the  case  of  superlattices,  the  overlap  of  wavefunctions  causes  a  splitting  of 
the  degenerate  energy  levels  in  the  neighboring  quantum  wells.  This  splitting  of 
discrete  energy  levels  results  in  the  creation  of  minibands.  The  minibands  consist  of 
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Figure  3  Energy  levels  and  minibands  in  a  multiple  quantum  well  structure  and 
a  superlattice  (19). 

a  number  of  individual  bands  equal  to  the  number  of  wells  in  the  superlattice  (15). 
The  presence  of  minibands  allows  a  variety  of  transitions  to  occur  corresponding 
to  different  wavelengths. 

In  the  case  of  multiple  quantum  well  structures,  the  lack  of  wavefunction 
overlap  leaves  the  energy  levels  discrete.  This  also  allows  an  increase  in  the  optical 
absorption  of  the  device  by  a  factor  equal  to  the  number  of  wells  in  the  structure. 
However,  ‘effective  minibands ’will  also  form  in  a  multiple  quantum  well  structure 
due  to  slight  variations  in  well  widths  or  alloy  composition  (15).  This  results  in 
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slightly  different  energy  levels  in  each  of  the  quantum  wells  making  up  the  multiple 
quantum  well  structure. 


1.2  Optical  Absorption  Properties 


The  foundation  of  optoelectronic  devices  lies  in  the  interaction  between  in¬ 
cident  electromagnetic  radiation  and  electrons  or  holes  in  the  detector  material. 
This  interaction  can  be  described  by  the  absorption  coefficient  a  (a  <  0  for  gain 
material  and  a  >  0  for  attenuating  material).  Absorption  occurs  when  a  photon 
assisted  transition  occurs;  an  electron  gaining  energy  from  a  photon  and  moving 
from  one  energy  level  to  another. 

The  form  of  the  linear  absorption  coefficient  used  by  Green  (1)  can  be  ex¬ 
pressed  as 


a  oc 


/  {Pnm  (k||)) 

En  (k||)  -  Em  (k,|) 


(1) 


where  Pmn  is  the  momentum  matrix  elements  for  transitions  between  states  n 
and  m  and  En  and  Em  correspond  to  the  energy  levels  of  the  initial  and  final 
states.  Both  the  energy  and  the  momentum  vector  are  functions  of  ky  and  must 
be  evaluated  at  each  value  of  ky  for  which  the  integral  must  be  performed.  The 
integral  in  equation  1  is  performed  through  the  region  of  k  -space  dictated  by  the 
Fermi  level. 


A  common  way  to  calculate  absorption  coefficients  is  to  perform  a  calcula¬ 
tion  of  the  quantum  well  band  structure  and  momentum  matrix  elements  for  one 
direction  in  k- space.  In  many  cases,  the  assumption  that  the  band  structure  is 
isotropic  throughout  the  plane  is  made  and  the  simple  factor  of  2tt  is  utilized. 
However,  the  band  structure  is  not  isotropic  throughout  the  plane  of  the  quantum 
well.  Therefore,  to  accurately  model  absorption  coefficients  the  magnitude  of  the 
anisotropy  must  be  determined  and  considered  in  any  calculation. 
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1.3  Anisotropy  Analysis 

There  are  several  models  used  for  calculating  the  bandstructures  for  semi¬ 
conductor  materials.  The  most  common  are  the  nearly  free  electron  model,  the 
tight-binding  method  and  the  one  used  in  this  analysis,  the  k  •  p  method  (17). 
The  nearly-free-electron  model  assumes  free  electrons  interacting  weakly  with  the 
positive  ions  in  the  lattice.  Due  to  Bragg  reflection  forbidden  energy  zones  occur, 
i.e.  bandgaps  develop.  In  the  tight-binding  model,  electronic  wave  functions  are 
assumed  to  be  comprised  of  hydrogenic  wavefunctions  which  are  localized  at  the 
lattice  sites.  The  potential  is  described  by  Coulombic  potentials  of  the  nearest 
neighbors  only.  Though  the  tight-binding  model  is  considered  the  most  accurate, 
it  is  computationally  intensive.  The  final  method  uses  the  k  •  p  model  which  will 
be  described  in  detail  in  chapter  2.  This  is  a  semi-empirical  method  based  on 
experimentally  measured  quantities. 

1.4  Previous  Work 

Over  the  past  five  years  modeling  of  the  absorption  spectra  of  quantum  well 
structures  has  been  the  subject  of  two  graduate  research  projects  at  the  Air  Force 
Institute  of  Technology  (AFIT)  (2).  Mathematical  models  have  been  developed 
and  experimental  data  has  been  taken  for  Si/SiGe  quantum  well  structures.  These 
models  allow  calculation  of  energy  levels  in  the  quantum  well  and  the  absorption 
coefficients  corresponding  to  bound-to-bound  valence  band  transitions.  The  cal¬ 
culations  were  based  on  the  band  structure  corresponding  to  a  single  azimuthal 
direction  in  the  plane  of  the  quantum  well  and  the  assumption  of  band  structure 
isotropy. 

As  discussed  earlier,  the  calculation  of  the  absorption  coefficients  requires 
an  integration  of  the  energy  bands  throughout  the  plane  of  the  quantum  well. 
This  requires  a  calculation  of  the  band  structure  for  each  incremental  step  in  the 
azimuthal  angle  theta.  This  research  reproduces  the  mathematical  models  of  the 
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quantum  well  structure  developed  by  Gregg  (2)  in  FORTRAN  code  and  allows 
enhanced  analysis  of  the  band  structure  and  the  associated  anisotropy. 

The  remainder  of  this  document  presents  the  theory  and  the  computational 
approach  used  to  determine  the  band  structure  in  a  semiconductor  quantum  well. 
First,  bulk  material  band  structures  are  calculated  using  the  Luttinger-Kohn  k  •  p 
theory  including  spin-orbit  interaction  effects.  The  strain  Hamiltonian  developed 
by  Kleiner  and  Roth,  presented  later  by  Hasegawa  (4)  and  then  by  People  and 
Jackson  (8)  is  also  incorporated.  Second,  the  analysis  for  the  quantum  well  is  pre¬ 
sented.  Quantum  well  valence  band  offsets  are  based  on  model  solid  theory  of  Van 
De  Walle  as  presented  by  People  and  Jackson  (8).  The  Envelope  Function  approxi¬ 
mation  as  presented  by  Szmulowicz  (6)  is  used  to  determine  the  quantum  well  band 
structure.  The  FORTRAN  programs  necessary  to  perform  the  required  calcula¬ 
tions  are  also  summarized  and  discussed.  Finally  the  resulting  dispersion  curves 
for  the  quantum  well  are  presented  for  various  azimuthal  angles  through  the  plane 
of  the  quantum  well.  Materials  considered  are  Si/Sii-zGe*  ,  GaAs/ Al,:Gai_xAs  , 
and  GaAs/InxGai-zAs. 
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II.  Theoretical  Approach 

This  chapter  discusses  in  general  the  approach  taken  in  determining  the  band  struc¬ 
ture  for  a  semiconductor  quantum  well.  Section  2.1  describes  the  development  of 
the  total  Hamiltonian  for  the  bulk  material  and  how  bulk  material  bandstructures 
are  determined  form  the  eigenvalues  of  the  Hamiltonian  matrix.  Section  2.2  dis¬ 
cusses  the  use  of  the  Envelope  Function  Approximation  (EFA)  in  the  modeling  of 
the  semiconductor  quantum  well  and  compares  its  development  and  solution  to 
that  of  the  one- dimensional  square  well  problem  solved  in  basic  quantum  mechan¬ 
ics. 

2. 1  Bulk  Material 

The  determination  of  energy  levels  in  a  bulk  material  semiconductor  is  ac¬ 
complished  by  solving  the  time-independent  Schrodinger  equation. 

Hip  =  Exp  (2) 

The  general  form  for  the  bulk  material  Hamiltonian  is  comprised  of  several 
terms.  These  are  the  k  •  p  Hamiltonian  Hkp,  the  spin-orbit  Hamiltonian  Hso,  and 
the  strain  Hamiltonian  H strain- 

2.1.1  k  ■  p  Theory.  Electrons  and  holes  in  a  periodic  potential  such  as  a 
crystal  lattice  obey  Bloch’s  Theorem.  Bloch’s  theorem  states  that  wavefunctions 
in  a  periodic  potential  can  be  expressed  as  the  product  of  a  function  having  the 
periodicity  of  the  lattice  and  a  plane  wave. 

The  general  form  for  such  a  function  can  be  expressed  as 

^nk  (r)  =  eik'ru„k  (r)  (3) 
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where  u  is  a  function  with  the  periodicity  of  the  lattice.  Substituting  this  into  the 
one-electron  Schrodinger  equation 


V'nk  =  En  (k)  V’nk 


(4) 


and  assuming  a  periodic  potential  will  give  the  general  k  •  p  equation(7) 


V  +  fik-P 

2m  m 


+ 


h2k2 
2m  ^ 


V(r) 


unk  (r)  =  En  (k)  unk  (r) 


(5) 


Bloch  functions  of  the  form  in  equation  3  create  a  complete  set  of  basis 
functions  for  any  value  of  k.  One  possibility  used  by  Luttinger  and  Kohn  is  based 
on  Bloch  functions  at  zone  center,  k  =  0  (3).  Following  an  approach  outlined  by 
Kane  (18),  un k  (r)  in  equation  5  can  be  expressed  in  the  k  •  p  representation  by 


Wnk  (r)  —  ^  ^  C-ri’n  Uji'ko  (r)  (®) 

n' 

Multiplying  both  sides  of  equation  5  by  unko  (r)  ,  replacing  unk  (r)  with  the  form 
given  in  equation  6,  and  integrating  gives 

&nn'  H  k  ■  p n'n  f  Cn'n  =  En  (k)  Cnn  (7) 

m  J 

The  term  pn/n  is  given  by 


En  (ko)  + 


fi2fc2 

2m 


P  n'n  — 


(r)  pu„/k0  (r)  dr 


(8) 


where  the  integral  is  over  the  unit  cell  and  p  is  the  momentum  operator  given  by 


p  =  —ihV 


(9) 
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Using  an  application  of  perturbation  theory  for  small  k  to  solve  equation  7 , 
Luttinger  and  Kohn  derived  the  general  form  for  the  Hamiltonian  describing  elec¬ 
tron  or  hole  behavior  in  a  periodic  potential  (3).  However,  the  momentum  matrix 
elements  given  in  equation  8  are  difficult  to  calculate  because  they  require  specific 
knowledge  of  the  Bloch  functions.  A  more  manageable  expression  of  the  Hamil¬ 
tonian  uses  angular  momentum  operators  and  inverse  mass  parameters  known  as 
Luttinger  coefficients  (2).  By  using  these  experimentally  measured  values,  the  cal¬ 
culation  of  the  momentum  matrix  elements,  p n>n,  of  equation  10  are  avoided.  Using 
these  values,  the  Hamiltonian  can  be  expressed  as 

tfk  p  =  2rn,  ^ xx  "1"  kZz)  (10) 

— 1272  —  kxx  +  c.p^j  —  24q3  (^(^LXXJ  kxy  +  c.p.^j 


where  q's  are  the  empirically  determined  Luttinger  coefficients  and  c.p.  is  a  cyclic 
permutation  in  x ,  y,  and  z.  Equation  10  gives  three  doubly  degenerate  energy 
bands  in  the  dispersion  relationship  which  are  degenerate  at  zone  center.  The  three 
bands  are  labeled  heavy  hole,  light  hole,  and  split  off  based  on  their  corresponding 
effective  masses. 

2.1.2  Coordinate  System  Transformation.  Since  the  previous  equations 
have  been  developed  using  symmetry  based  coordinates,  a  transformation  must  be 
done  to  align  the  coordinate  system  with  the  direction  of  interest.  This  is  done 
through  the  use  of  a  rotation  matrix  based  on  euler  angles  applied  to  both  the 
angular  momentum  operators  and  the  wave  vectors  (1). 
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The  transformation  is  given  as 


where  R  is  given  by 


cos  9  cos  <f>  —  sin  <j>  sin  9  cos  (ft 
R  =  cos  9  sin  4>  cos  </>  sin  9  sin  (f> 
i  —  sin#  0  cos# 


(12) 


(13) 


The  values  for  9  and  <j>  are  based  upon  the  direction  of  crystal  growth.  For  the 
[110]  direction  9  =  §  and  <j>  =  f .  For  the  [001]  direction  9  =  0  and  </>  = 

2.1.3  Spin- Orbit  Hamiltonian.  Diagonalizing  the  form  of  the  Hamiltonian 
given  in  equation  10  results  in  a  six  fold  degeneracy  at  zone  center.  If  the  interaction 
between  the  electron’s  spin  and  orbital  angular  momentum  is  incorporated,  the 
degeneracy  will  be  partially  removed.  Terms  describing  the  spin-orbit  induced 
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perturbations  are  diagonal  elements  in  the  \j,  rrij,  k)  basis  and  are  given  by 

0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  -A  0 

0  0  0  0  0  -A 

where  A  is  the  experimentally  measured  spin-orbit  interaction  term,  j  = 

™  -  3  1  _1  _3  flTlr|  a  _  I  m  .  _  1  _1 

m3  •  2 5  2  ’  2 5  2  dI1U  3  ~  2’  '2’  2’ 

k  ■  p  Hamiltonian.  Performing  the  transformation  described  in 
section  2.1.2  results  in  the  following  form  for  the  bulk  material  k  •  p  Hamiltonian 
in  the  \j,mj,k)  basis: 

hh  a  0  0  ij=  -ipV 2 

a*  Ih  0  0  iajl 

0*  0  Ih  -a  *75 

0  0*  -a*  hh  -i0*V 2 

ic*\/l  ^  50  0 

0-V2  -i%  0  so 
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For  the  [001]  growth  direction: 


hh  =  (kl  +  kl  +  kl)y1  +  (kl  +  kl-2klh2  (16a) 

Ih  =  (£  +  +  fcfrh -(*?  +  *$- 2*Jbs  (16b) 

so  =  (kl  +  kl  +  kl  )Tl  (16c) 

a  =  —2i'JZ{kx  —  iky)kznfi  (16d) 

/3  =  -  k’’-}'/,  -  2ikJ;,~i:t  (16e) 

6  =  Ih-hh  (16f) 


For  the  [110]  growth  direction: 


hh  = 

{k2x  +  ky  +  k2z)  7i  —  {k2  —  2k2  +  +  (A '2  —  k2)^  73 

(16g) 

Ih  = 

(kl  +  k2  +  k2)  +  ( k2  —  2k2  +  k2z )y  —  ( k \  —  k2)-  73 

(16h) 

so 

C kl  +  k2  +  k2z) 

(16i) 

<\  = 

V^3(^aj/y2  ^/^y/y 

(16j) 

(3  = 

y  (( fc*  -  +  /cf)72  +  -  tikxky  -  &z)73) 

(161c) 

8  = 

//i  —  /l/l 

(161) 

2.1.5  Strain  Hamiltonian.  When  two  materials  having  different  lattice 
constants  are  used  to  create  a  heterostructure  ,  strain  will  be  induced  into  the 
material  as  a  result  of  the  necessary  compensation  for  the  lattice  mismatch.  This 
perturbation  to  the  symmetry  of  the  periodic  potential  changes  the  energy  levels 
in  such  a  way  as  to  cause  a  splitting  in  the  degeneracy  of  the  light  hole  and  heavy 
hole  bands.  The  Hamiltonian  describing  the  effect  of  strain  is  given  in  terms  of 
orbital  angular  momentum  operators;  deformation  potentials,  Dy ,  D2,  and  D3;  and 


strain  tensor  components  eap  (2) 


H strain  —  (®ix  d"  &yy  d"  &zz)  I  (1^) 

+2D2  ^ L2X  -  exx  +  c.p.  +  2D3  (t^  exy  +  c.p. 

where  c.p.  denotes  cyclic  permutation.  The  first  term  of  equation  17  represents  a 
shift  in  average  energy  of  the  valence  band.  Therefore,  this  term  is  excluded  from 
bulk  calculations  since  absolute  values  of  energy  are  not  an  issue.  However,  in  the 
quantum  well  band  structure  calculations,  where  absolute  or  relative  energy  values 
are  important,  this  term  is  incorporated. 

Performing  a  rotation  operation  on  the  angular  momentum  operators  of  equa¬ 
tion  17  together  with  a  similarity  transformation  of  the  strain  tensor  similar  to  that 
done  earlier  for  Hk  p,  yields  the  strain  Hamiltonian  where  the  Dx  term  in  equation 
17  has  been  omitted 


e2 

0 

e3 

0 

0 

—ie3\J 2 

0 

— e2 

0 

e3 

-ie2\/ 2 

0 

e3 

0 

-e2 

0 

0 

-ie2V 2 

0 

e3 

0 

e2 

-ie3V 2 

0 

0 

ie2V 2 

0 

ie3\f 2 

0 

0 

ie3V 2 

0 

ie2yj 2 

0 

0 

0 

(18) 


For  the  [001]  growth  direction: 

(19a) 
(19b) 


^2  —  qD2(Czz  &xx) 

e3  —  0 
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Figure  4  Compensation  for  lattice  mismatch  in  a  strained  material  (2)  . 

For  the  [110]  growth  direction: 

62  =  ^  £*2  (fizz  &xx)  d-  "^D^eXy  (20a) 

1  1 

63  =  — ^=Z?2(6xx  g22)  +  -^^=D^6Xy  (20b) 

The  values  for  D{  are  based  on  linear  interpolation  of  the  deformation  po¬ 
tentials  for  the  alloy’s  constituents.  D\  is  the  interpolated  value  of  the  hydrostatic 
deformation  potentials  av,  £>2  is  the  interpolated  value  of  Du,  the  valence-band 
deformation  potential  for  uniaxial  stress  in  the  [001]  direction,  and  £>3  is  the  inter¬ 
polated  value  of  D'u,  the  valence  band  deformation  potential  for  uniaxial  stress  in 
the  [111]  direction. 

The  values  eap  are  conventional  strain  components  and  can  be  expressed  in 
terms  of  the  unique  elements  of  the  stiffness  matrix  C  for  cubic  structures  and  the 
value  for  in-plane  strain  [(8), (7)].  The  in-plane  strain  is  a  function  of  the  lattice 
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constants  in  the  well  and  the  barrier  material  and  can  be  expressed  as 


(21) 


where  b  is  the  unstrained  lattice  constant  and  as  is  the  lattice  constant  of  the 
substrate  or  barrier  material.  This  represents  the  percentage  change  in  the  lattices 
constant  (7). 

The  expressions  for  the  strain  components  are  (8) 
for  the  [001]  growth  direction: 


&XX 

&xy 

&zz 


eu 

0 


for  the  [110]  growth  direction: 

&XX  =  &\\ 


'LJ  [2^44  +  Cn  +  C\2 

&zz  =  &\\ 


2C44  —  Cvi 
2O44  +  Cn  +  C\2 
2  [Cn  —  2C12] 


(22) 


(23) 


2.1.6  Energy  Level  Band  Structure  for  Bulk  Material.  The  Hamiltonian 
for  unstrained  bulk  material  can  now  be  written  as  Htotai  =  Hkp  +  Hso  and  the 
Hamiltonian  for  strained  bulk  material  can  be  written  as  Htotai  =  Hk-P  +  Hso  + 
H strain-  Both  expressions  are  functions  of  k  =  (kx,  ky,  kz)  =  (/c*,  ky)  and  empirically 
derived  or  tabulated  data  specific  to  the  constituent  elements.  These  constants 
are  listed  in  appendix  B.  Using  the  approach  by  Gregg  (2)  ,  the  values  used  in 
the  calculation  for  an  alloy  material  are  linearly  interpolated  values  based  on  the 
percentage  composition  of  one  material  in  the  other.  For  example,  the  Luttinger 
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coefficients  in  an  alloy  of  Silicon  and  Germanium  are  given  as 

7  =  7si  +  *(7 Ge  -  7si)  (24) 

where  x  is  the  percentage  of  Germanium  in  the  alloy. 

Diagonalizing  the  appropriate  Hamiltonian  gives  the  energy  values  for  the 
bulk  material  as  a  function  of  k.  In  an  unstrained  material,  there  will  be  two  bands 
at  zone  center,  one  being  doubly  degenerate,  the  other  being  quadruply  degenerate. 
The  doubly  degenerate  band  is  the  split  off  band  and  the  quadruply  degenerate 
bands  at  zone  center  are  the  light  hole  and  heavy  hole  bands.  The  bands  are 
named  for  their  effective  masses  and  their  relative  location  as  a  function  of  energy. 
In  strained  material  there  will  be  three  doubly  degenerate  bands.  This  is  a  result 
of  the  light  hole  and  heavy  hole  bands  splitting  due  to  the  strain  perturbation. 
Examples  of  the  bandstructure  are  shown  in  figure  5  for  Si/Si7oGe3o  alloy. 


2.2  Quantum  Well 

The  approach  used  in  determining  the  energy  eigenstates  within  a  semicon¬ 
ductor  quantum  well  structure  is  similar  to  that  used  in  determining  the  eigenstates 
for  a  one  dimensional  square  well  potential.  This  section  first  revisits  the  one  di¬ 
mensional  square  well  potential  and  then  explains,  in  an  analogous  manner,  the 
solution  to  the  semiconductor  quantum  well  problem. 


2.2.1  One  Dimensional  Square  Well  Problem  (20).  Recalling  the  sim¬ 
ple  one  dimensional  square  well  problem,  one  starts  with  the  time  independent 
Schrodinger  equation 


+  V(z) 


xj;  —  E'lp 


(25) 
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(a)  (b) 

Figure  5  Bulk  material  dispersion  relations  for  Si/Si7oGe30  resulting  from  diag- 
onalization  of:  (a)  Hfc  p  +HS0  and  (b)  Hfcp  +  Hso  +  Hs(roin 


where  p  is  the  momentum  operator  given  as 

,Kd 


(26) 


Next  the  wave  functions  are  assumed  to  be  linear  combinations  of  plane  waves  of 
the  form 


(Z)left 

=  AekbZ  +  Be~kbZ 

(27a) 

(Z)well 

_  Qgikwz  _|_  JJg~ikwz 

(27b) 

^  (^0  right 

=  EekbZ  +  Fe~kbZ 

(27c) 

Using  this  form  for  the  wave  functions  in  the  Schrodinger  equation  given  in  equation 
25,  values  for  the  wave  vector  k  can  be  expressed  as  a  function  of  energy  E.  The 
expressions  for  k  in  the  well  and  barrier  of  a  simple  square  well  potential  are 


ku 


h 


\/2 rnE 

h 

y/2m(Vr0  -  E) 


h 


(28a) 

(28b) 


Since  the  wave  functions  are  utilized  to  describe  the  behavior  of  observables  such 
as  position,  momentum,  etc.,  boundary  conditions  must  be  placed  on  the  wave 
functions  at  the  well/barrier  interface  to  ensure  both  the  wave  functions  and  their 
derivatives  are  well  behaved.  These  boundary  conditions  are 


i’(4)ri’(4)w  (29a) 

JjV'WL*!  =  |VmL±.  (29b) 

An  addition  boundary  condition  placed  on  the  wave  functions  is  that  they  must 
decay  to  zero  as  z  approaches  +oo.  This  requires  the  coefficients  B  and  E  to  be 
equal  to  zero.  Using  these  boundary  conditions  at  the  barrier/well  interface,  a 
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system  of  equations  can  be  developed  and  expressed  as 


^e-ifcM(f)  _|_  £)eik™(  2) 

=  Ae  ^(2) 

(30a) 

Ceikw(  a)  +  De~ikw(*) 

=  Fe~kb& 

(30b) 

ikwCe~ikw(%)  -  ikwDeikw& 

=  kbAe~h (*) 

(30c) 

ikwCeikw (t)  -  ikwDe~ikw (2) 

=  - kbFe~kb (t) 

(30d) 

This  system  of  equations  can  be  represented  by  a  matrix  equation  of  the  following 
form 

e-ifcu,(|)  e-ifc“(f)  _e~M§)  o  \  /  \  0 

e~ikw{%)  e-ife«,(f)  Q  _e-fci>(f)  £)  Q 

ikwe~ikw (§)  -ikwe~ikw(*)  - kbe~kb (a)  0  .4  0 

-ikwe~Mi)  0  fc6e~fc6(t)  )  \  F  )  \0 

or  more  simply 

[M]  -  x  =  0  (32) 

This  can  only  be  true,  assuming  non-trivial  elements  of  Xi  if 

det  [M]  =  0  (33) 

The  matrix  M  is  a  function  of  k  which  is  a  function  of  energy  E.  Thus 
an  expression  for  det  [M]  as  a  function  of  energy  can  be  derived  and  the  roots 
found  with  a  numerical  root  finding  tool.  These  roots  correspond  to  the  energy 
eigenvalues  of  the  system. 

2.2.2  Semiconductor  Quantum  Well.  The  reason  for  the  simplicity  of 
the  calculation  in  the  simple  square  well  problem  is  that  k  can  be  found  easily 
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Figure  6  One  dimensional  simple  square  well  potential  of  width  a. 

as  a  function  of  energy.  For  the  case  of  semiconductor  quantum  well  structures 
this  relationship  cannot  be  expressed  analytically.  Therefore  a  method  must  be 
developed  to  invert  the  expressions  found  in  the  previous  chapter  from  E  (k)  to 
k  ( E )  or  more  precisely  kz(E ,  ky).  This  is  accomplished  through  a  method  presented 
by  Szmulowicz  which  expands  the  Hamiltonian  in  powers  of  kz  and  develops  an 
eigenvalue  problem  from  which  kz  can  be  found  for  a  particular  set  of  values  for 
energy  E  and  a  particular  value  for  ky  (6).  This  process  is  explained  in  section 
2. 2.2. 2. 

Once  the  values  for  kz  are  known,  they  can  be  incorporated  into  a  matrix 
equation  of  the  form  shown  in  equation  33  to  see  if  the  determinant  is  equal  to  zero. 
Recall  that  the  matrix  was  developed  from  the  application  of  boundary  conditions 
to  the  well  and  barrier  wave  functions.  If  the  determinant  is  equal  to  zero,  then  the 
values  of  E  are  in  fact  eigenvalues  for  the  system.  Using  this  approach,  dispersion 
relationships  can  be  developed  for  the  quantum  well  structure.  This  process  is 
explained  in  sections  2. 2. 2. 3  through  2. 2. 2. 4. 

2.2.2. 1  Valence  Band  Energy  Offset.  For  the  bulk  material  calcu¬ 
lations  in  the  previous  section  of  this  chapter,  absolute  values  for  energy  were  not 
considered.  The  only  objective  was  to  present  the  bands  relative  to  one  another. 
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However,  for  the  calculation  of  the  quantum  well  it  is  necessary  to  predict  not  only 
an  absolute  value  for  the  bands  but  also  the  difference  or  offset  between  bands  in 
the  barrier  and  the  well  material.  This  is  analogous  to  the  depth  of  the  potential 
well  in  the  one  dimensional  problem. 

The  approach  taken  for  this  model  uses  an  average  valance  band  energy  E”v 
based  on  model  solid  theory  of  Van  De  Walle  and  presented  by  People  and  Jackson 
(8).  Model-solid  theory  was  developed  to  estimate  the  electrostatic  potential  differ¬ 
ence  across  a  heterojunction.  The  value  of  E%v  corresponds  to  the  average  energy 
at  the  top  of  the  valence  band,  ignoring  spin-orbit  interaction  effects.  Using  this 
basis  to  establish  an  absolute  scale,  the  average  energy  for  the  valence  band  in  a 
material  is  a  function  of  three  values:  the  average  valence  band  energy  Ehav  found 
through  model  solid  theory;  the  magnitude  of  the  spin-orbit  term  A;  and  the  shift 
in  the  average  energy  due  to  crystal  strain. 

The  first  two  values  can  be  found  tabulated  (8).  The  shift  due  to  the  strain 
is  calculated  from  tabulated  data  and  is  a  function  of  an  alloy’s  stoichiometry.  The 
shift  in  the  average  valence  band  energy  due  to  strain  is 

=  D\(2exx  +  ezz)  (34) 

where  D\  is  a  constant  based  on  linear  interpolation  of  the  hydrostatic  deformation 
potentials  av  for  the  two  elements  in  the  alloy.  The  expressions  for  exx  and  ezz  were 
given  earlier  in  equations  22  and  23. 

Mathematically,  the  offset  is  incorporated  by  adding  a  constant  to  both  the 
well  Hamiltonian  and  the  barrier  Hamiltonian.  Therefore  the  total  Hamiltonians 
are  expressed  as 


H total  —  Hk-p  "b  Hso  +  Vb  Eshift 

H total  =  Hkp  ~\~  H so  -\~  H strain  H""  Uu  E shift 
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(35a) 

(35b) 


The  values  for  Vi ,  and  Vw  are  determined  based  on  the  average  valence  band 
energies  of  the  well  and  barrier  material.  In  order  to  stay  consistent  with  conven¬ 
tion,  a  constant  value  is  included  that  will  shift  the  zone  center  (k  =  0)  heavy  hole 
band  in  the  well  bulk  material  to  zero.  The  values  for  Vj,  and  Vw  can  be  expressed 
as 


Vw  = 
Vb  = 
Eshift  = 


A  well 

tp  well  ,  ^  |  _ well 

Kyav  +  -3“  +  el 


rpbarr  , 
^v.av  1 


A  well 

TP  well  1  „well  1  ^  1  tp 

K,av  +  ei  +  — r~  +  El 


(36a) 

(36b) 

(36c) 


The  term  Ehh  comes  from  the  diagonalization  of  the  matrix  Hso  +  Hstrain 
and  represents  the  shift  in  the  heavy  hole  band  from  the  unstrained  position  in  the 
bulk  dispersion  relations.  Again,  the  objective  of  subtracting  E shift  is  to  bring  the 
zone  center  energy  of  the  heavy  hole  band  in  the  bulk  dispersion  relations  of  the 
well  material  to  zero  and  is  done  for  convention  only. 
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2.2. 2.2  kz  as  a  Function  of  Energy  and  k-parallel.  Now  that  the 
Hamiltonians  for  both  the  well  and  the  barrier  are  well  defined  for  any  value  of  k, 
the  next  step  is  to  invert  the  expression  to  a  from  of  kz  as  a  function  of  energy  E 
and  k||  .  This  is  done  through  a  method  presented  by  Szmulowicz  (6). 

The  total  Hamiltonian,  given  by  Hk.p  +  Hs0  +  Hstrain  can  be  factored  in  terms 
of  kz  where 


H  =  H^k^  +  H\kz  +  Hq 


(37) 


and  #2,  Hi,  and  H0  are  6x6  matrices  that  are  functions  of  ky  =kx  +  ky.  Incorpo¬ 
rating  this  expression  of  the  Hamiltonian  into  the  time-independent  Schrodingei 
equation 

Hip  =  E4 >  (38) 

gives  the  following  expression  relating  kz,  E,  H2,  Hi,  and  H0 

k\4)  +  H?Hxkz4>  +  Hz1  (H0-E)4>  =  0  (39) 


Defining  a  vector  (f>  as 


(j)  =  kzif 


(40) 


and  substitution  into  equation  39  yields 


kz4>  +  H^HX(t>  +  tf2_1  {Ho-E)ip  =  0  (41a) 


kz4>  -  (j>  =  0 


(41b) 
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Through  algebraic  manipulation,  the  system  of  equations  can  be  expressed  in  a 
matrix  equation  of  the  following  form 


(  -kzI  I 

\  - H 71  (H0  -  E)  -H?HX  -  kzI 


(42) 


where  I  is  a  6x6  unit  matrix.  If  kz  is  considered  an  eigenvalue,  the  problem  can  be 
considered  an  eigenvalue  problem  of  the  form  (M  —  XI)  x  —  0  expressed  as 


Y  0 

\-H^(H0-E)  -H?HX 
where 


kzI  0 
0  kzIex6 


M  = 


0  Isx6 

-H.;1  (H„  -  E) 


(44) 


and 

By  finding  the  eigenvalues  of  this  complex  non-Hermitian  matrix,  kz  can  be  ex¬ 
tracted  for  any  value  of  energy  E.  By  finding  the  eigenvectors  of  the  matrix  M,  the 
eigenvectors  associated  with  each  value  of  kz  can  also  be  determined.  The  method 
used  to  determine  the  eigenvectors  and  eigenvalues  are  presented  in  section  3.2. 


2. 2. 2. 3  Envelope  Functions.  Envelope  functions  are  linear  com¬ 
binations  of  bulk  material  wavefunctions  where  the  periodic  piece  of  the  Bloch 
function,  un k=o  (r) ,  is  not  included.  The  flat-band  approximation  states  that  the 
periodic  part  of  the  Bloch  functions  in  the  well  and  the  barrier  are  approximately 
equal,  i.e.  un k=o  (r)  is  continuous  (13).  Using  this  approximation,  «nk=o  (r)  can  be 
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excluded  when  imposing  boundary  conditions.  Through  symmetry  analysis,  the 
envelope  functions  in  the  well  region  and  the  barrier  region  can  be  expressed  in 
matrix  notation  as 


Fb  (k||,  z)  =  (TCb)Eb(-z)B  (46) 

Fw{khz)  =  [CwEw{z)±TCwEw(-z)]W  (47) 


where  T  is  a  reflection  operator  defined  as 

^  1  0  0  0  0  0  ^ 

0-100  00 

0  0  1  0  0  0 

r  = 

0  0  0  -1  0  0 

0  0  0  0  -1  0 

^  0  0  0  0  0  1  ) 


(48) 


and  represents  a  reflection  about  the  origin.  This  operator  is  used  with  the  parity 
of  the  wave  functions  to  reduce  the  number  of  boundary  conditions  (2). 

C  is  a  matrix  consisting  of  the  eigenvectors  specifying  the  linear  combination 
of  plane  waves  making  up  the  bulk  material  wave  functions.  The  matrix  C  consists 
of  the  eigenvectors  in  column  form  and  is  determined  by  the  process  described  in 
the  previous  section.  C  can  be  expressed  as 


c  =  (  i>hh  i’lh  4>ih  hh  $ so  4>. so )  (49) 

where  ipx  represent  the  eigenvectors  associated  with  the  eigenvalues  of  kz  in  column 
form. 

The  E  matrices  contain  the  plane  wave  portions  of  the  wave  functions  and 
utilize  the  eigenvalues  of  kz.  They  are  defined  in  both  directions.  For  the  positive 
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z  direction 


eik‘ *hz  0  0  0  0  0 

0  eik*hz  0  0  0  0 

0  0  eik*hz  0  0  0 

0  0  0  eik*hz  0  0 

0  0  0  0  eifc*°2  0 

0  0  0  0  0  .  eik‘°z 

and  for  the  negative  2  direction 

e-ik*hz  0  0  0  0  0 

0  e-iklzhz  0  0  0  0 

0  0  e-ik‘zhz  000 

0  0  0  e~ik*hz  0  0 

0  0  0  0  e~ik°°z  0 

0  0  0  0  0  e~ik>°z 

W  and  B  are  eigenvectors  that  specify  the  linear  combination  of  bulk  ma¬ 
terial  wavefunctions  that  make  up  the  envelope  functions  in  the  well  and  barrier 
respectively. 

2.2.24  Boundary  Conditions.  The  boundary  conditions  are  similar 
to  those  required  in  the  one-dimensional  square  well  problem.  For  the  simple 
square  well  problem  the  boundary  conditions  require:  (1)  the  wave  functions  are 
continuous  at  the  well/barrier  boundaries,  (2)  the  derivatives  of  the  wave  functions 
are  continuous  at  the  well/barrier  boundaries,  and  (3)  the  wave  functions  approach 
zero  at  +00.  These  wave  functions  and  their  derivatives  must  be  continuous  in 
order  to  correspond  to  an  observable  probability  density.  If  the  derivatives  were 
not  continuous  the  presence  of  a  singular  potential  would  be  implied.  Finally,  the 


wave  functions  must  be  finite  everywhere  and  square  integrable  which  requires  the 
wave  functions  to  approach  zero  at  +oo  (14). 

The  requirements  for  the  semiconductor  quantum  well  are  essentially  the 
same  as  those  stated  above.  For  bound  carriers,  the  wavefunctions  must  decay 
exponentially  in  the  barrier  to  ensure  the  wave  functions  remain  finite  everywhere. 
The  wave  functions  must  also  be  continuous  at  the  heterostructures  as  before. 
However,  due  to  the  differences  in  bandstructure  in  the  barrier  and  the  well  ma¬ 
terials,  the  effective  masses  will  change  at  the  barrier/well  interface.  For  a  simple 
ID  quantum  well  this  results  in  the  boundary  condition  for  the  derivative  to  be 

—  —Fb=——Fw  (52) 

mb  dz  mw  dz 


Through  the  use  of  the  definition  of  the  effective  mass  and  the  expression  of  the 
total  Hamiltonian  factored  with  respect  to  kz,  the  boundary  conditions  using  the 
6x6  k  ■  p  bulk  eigenfunctions  can  be  expressed  as  (2) 


pb  __  pu 


rib  ®  pb  ,  pb 

HiTzF  +  - 2iF 


B  Hw 

Ho  zr~Fw  H - lr.Fw 

dz  -2i 


(53) 

(54) 


where  Fb  and  Fw  are  the  envelope  functions  in  the  barrier  and  the  well  and  Ho 
and  Hi  are  from  equation  37. 

Using  the  boundary  conditions  stated  in  equations  53  and  54  and  the  form  of 
the  envelope  functions  in  equations  46  and  47,  a  matrix  equation  can  be  expressed 
as 


Mw 

H?NW  +  ^Mu 


-Mb 
-( H%Nb  + 
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where  the  values  of  Mw ,  Mb,  Nw,  Nb  are  defined  as 


Mw 

=  CwEw{z)  +  TCwEw(-z) 

(56a) 

Mb 

=  CbEb{z ) 

(56b) 

Nw 

=  CWKWEW  {-z)  -  TCWKWEW  {. z ) 

(56c) 

Nb 

=  CbKbEb(z ) 

(56d) 

The  matrices  C,  T,  E  (- z ) ,  and  E  (+z)  were  defined  in  section  2. 2.2. 2 

The  K  matrices  consist  of  the  eigenvalues  for  kz  =  kz(E)  found  earlier  from 
the  eigenvalues  of  equation  44  and  are  arranged  as  follows  for  both  the  well  and 
barrier 

kbh  0  0  0  0  0 

0  klzh  0  0  0  0 

0  0^000 

0  0  0  kbh  0  0 

0  0  0  0  ksz°  0 

0  0  0  0  0  fcf 

where  hh,  Ih,  and  so  correspond  to  the  heavy  hole,  light  hole,  and  split  off  bands. 

2.2.3  Quantum  Well  Dispersion  Relationships.  As  was  done  for  the  sim¬ 
ple  quantum  well  problem,  the  matrix  equation  describing  the  boundary  conditions 
given  in  equation  55  can  be  assumed  to  be  an  eigenvalue  problem,  with  an  eigen¬ 
value  equal  to  zero.  Therefore,  the  boundary  conditions  can  only  be  satisfied  if  the 
determinant  of  the  matrix  in  equation  55  is  equal  to  zero.  The  methods  used  to 
compute  the  matrices  described  in  this  section  and  ultimately  find  the  determinant 
of  the  boundary  conditions  matrix  as  a  function  of  energy  E  are  described  in  the 
following  chapter. 
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III.  Computational  Approach 

This  chapter  explains  the  approach  taken  in  calculating  dispersion  curves  for  both 
bulk  material  and  quantum  wells.  Also  discussed  is  the  development  of  the  com¬ 
putational  tools  necessary  to  perform  the  specific  calculations. 

3.1  Bulk  Material  Calculations 

The  bulk  material  band  structure  requires  determining  the  eigenvalues  of  a 
6x6  complex  Hermitian  matrix.  The  total  Hamiltonian  matrix  consists  of  the  k  •  p 
Hamiltonian,  the  spin-orbit  Hamiltonian,  and  the  strain  Hamiltonian.  The  system 
must  be  solved  for  every  value  of  the  inverse  wave  vector  k.  This  is  accomplished 
through  the  use  of  standard  eigenvalue  subroutines  discussed  in  section  3.3. 

3.2  Quantum  Well  Calculations 

Developing  the  quantum  well  dispersion  curves  first  requires  determination 
of  the  zone  center  energies  for  each  band.  This  is  accomplished  through  thorough 
analysis  of  the  function  relating  the  determinant  of  the  boundary  condition  matrix 
to  energy.  Each  root  of  this  complex  function  identifies  a  quantum  well  energy 
band  at  a  particular  value  of  ky .  These  values  are  then  utilized  as  initial  guesses  in 
a  root-finding  algorithm  that  determines  roots  of  the  function  for  each  incremental 
step  in  k|| . 

The  two  significant  steps  involved  in  the  quantum  well  problem  are  the  expres¬ 
sion  of  kz  as  function  of  energy  E  and  then  finding  the  determinant  of  the  matrix 
describing  the  quantum  well  boundary  conditions.  Both  steps  are  explained  in  the 
following  sections. 

3.2.1  Eigenvalues  and  Eigenvectors  of  kz(E).  The  total  Hamiltonian  is 
first  presented  in  the  factored  form.  Recall  from  section  2. 2. 2.2  that  the  Hamil- 
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Figure  8  Schematic  of  a  bulk  dispersion  curve  showing  energy  as  a  function  of 
kz  and  kz  as  a  function  of  energy  (2). 

tonian  is  factored  in  powers  of  kz  to  yield  Htotai  =  H2kl+Hj kz  +  H0.  The  coefficients 
of  the  factored  Hamiltonian  are  then  utilized  to  develop  the  matrix  equation  de¬ 
scribing  the  inverse  wave  vector  component  kz  =  kz  (E,  ky)  as  a  function  of  energy 
and  k||.  The  complete  set  of  12  eigenvalues  and  eigenvectors  for  this  system  con¬ 
sist  of  only  6  linearly  independent  vectors  and  6  doubly  degenerate  eigenvalues. 
Therefore,  selection  of  the  correct  eigenvalues  and  the  corresponding  eigenvectors 
is  necessary. 

If  X  is  the  transformation  matrix  used  to  diagonalize  the  matrix  in  equation 
43,  it  consists  of  the  system’s  eigenvectors  in  column  form  and  can  be  expressed  as 

*-(*“*”)  («) 

V  X21  X22  ) 

The  column  vectors  in  X2i  and  X22  will  be  equal  to  the  X\\  and  X\2  multiplied 
by  the  associated  eigenvalue.  This  is  apparent  through  the  expression  given  in 
equation  42.  As  shown  in  figure  8,  for  every  positive  value  of  kz  there  exists 
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a  corresponding  negative  value.  The  eigenvectors  associated  with  the  negative 
eigenvalues  can  be  shown  to  be  equal  to  the  vectors  associated  with  the  positive 
eigenvalues  with  the  exception  of  a  phase  factor.  Thus,  all  positive  eigenvalues 
and  associated  eigenvectors  are  selected  and  the  rest  excluded  from  further  cal¬ 
culations.  It  must  also  be  noted  that  kz  will  become  pure  imaginary  when  the 
value  of  energy  is  greater  than  the  maximum  energy  of  the  associated  band.  In 
this  case,  there  is  still  a  negative  pure  imaginary  eigenvalue  associated  with  every 
positive  imaginary  value.  The  linear  dependence  on  the  eigenvectors  hold,  so  only 
the  positive  imaginary  values  and  associated  eigenvectors  are  utilized. 

A  subroutine  must  select  all  positive,  pure  real  and  pure  imaginary,  eigen¬ 
values  and  associated  eigenvectors.  It  is  also  important  that  for  each  iteration  the 
list  of  eigenvalues  and  eigenvectors  must  be  ordered  the  same  way  before  being- 
passed  to  the  following  subroutine  for  construction  of  the  quantum  well  boundary 
condition  matrix.  The  form  of  the  matrices 


eigenvalues  = 


and  for  the  matrix  of  eigenvectors  C 


is  given  below. 


t  hh  ^ 

Ih 
Ih 
hh 

SO 

\ so  / 


(59) 


c  =  ( 'Phi,  4’ih  i’ll,  i>hh  i>,c ) 


(60) 


where  ipx  represent  the  eigenvectors  in  column  form. 

In  the  simplest  case,  ky  =  0,  the  bulk  dispersion  curves  behave  as  described 
above  and  as  shown  in  figure  8.  In  this  case,  the  identification  of  hh,  Ih,  and  so 
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eigenvalues  is  easily  accomplished.  However,  as  ky  increases,  choosing  the  correct 
eigenvalues  becomes  more  complicated  because  the  imaginary  values  of  kz  are  much 
less  predictable.  Also,  as  ky  increases  the  bands  become  double  valued  as  a  function 
of  energy.  In  this  case,  an  algorithm  is  forced  to  choose  only  one  of  the  two 
valid  solutions  and  assign  the  second  value  erroneously  to  another  band.  This 
situation  occurs  whenever  any  one  of  the  bands  has  a  maximum  not  at  ky  =  0. 
The  difficulties  in  assigning  eigenvalues  to  a  specific  band  can  be  seen  in  figures  9 
thru  11  which  show  the  bandstructure  determined  form  the  kz(E)  inversion  process. 
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Figure  12  Plot  of  the  determinant  of  the  30 A  Si/SiGe  quantum  well  boundary 
condition  matrix  vs.  energy.  kx=ky=0. 

3.2.2  Determinant  of  Boundary  Condition  Matrix.  Once  the  appropriate 
eigenvalues  and  associated  eigenvectors  for  kz  are  determined  the  matrix  equation 
describing  the  boundary  conditions  for  the  quantum  well  is  generated.  As  shown 
in  section  2.2.3,  an  acceptable  value  of  energy  is  found  for  a  particular  value  of  ky 
when  the  determinant  of  the  matrix  in  equation  55  is  equal  to  zero. 

An  analysis  of  a  plot  of  the  real  and  imaginary  values  of  the  determinant  of 
the  aforementioned  matrix  will  allow  identification  of  the  roots  of  the  function  for 
a  particular  value  of  ky.  Examples  of  these  plots  are  shown  in  figures  12  and  13. 
Figures  12  shows  a  wide  range  of  values  for  energy  but  closer  analysis  shows  that 
roots  only  exist  in  the  small  region  corresponding  to  the  well.  An  enhanced  plot 
of  the  data  for  the  region  of  the  well  is  shown  in  figure  12. 
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Imaginary 


Figure  13  Plot  of  the  determinant  of  the  30A  Si/SiGe  quantum  well  boundary 
condition  matrix  vs.  energy.  kx=ky=0. 
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Using  these  roots  as  seed  values,  the  process  is  repeated  at  increasing  values 
of  k||  with  a  root  finder.  ky  is  incrementally  increased  in  a  particular  direction 
in  the  plane  of  the  quantum  well  to  determine  the  zeros  at  that  new  value  of  ky. 
These  values  are  then  used  as  the  new  seed  values  for  the  next  iteration. 

Since  the  initial  roots  at  ky  =  0  were  determined  with  good  accuracy  in 
the  first  step  of  the  problem  analysis,  a  modified  secant  method  is  acceptable 
when  finding  roots  at  higher  values  of  ky .  The  secant  method  is  a  modification  of 
Newton’s  method  and  is  shown  graphically  in  figure  14. 

The  basic  equation  governing  the  secant  method  is 


x2  =  xi  -  f  (xi) 


Xi  -  Xp 

f  (*i)  -  /  (*o) 


(61) 


Enhancements  were  made  to  the  basic  algorithm  to  handle  difficult  roots  such  as 
those  shown  in  figures  15  and  16.  In  many  cases  such  as  that  shown  in  figure  15 
the  real  part  of  the  complex  function  or  the  imaginary  part  of  the  function  can  be 
approximately  zero.  The  root  finder  must  determine  the  best  choice,  either  the  real 
part  or  the  imaginary  part  of  the  function,  to  utilize  in  the  calculation  given  in  the 
equation  above.  This  is  accomplished  by  comparing  the  values  of  the  determinant 
in  the  vicinity  of  the  seed  values.  The  absolute  values  of  the  real  and  imaginary 
part  are  compared.  If  one  is  approximately  zero  within  a  tolerance,  the  other  is 
used  for  the  calculation. 


In  other  cases,  such  as  that  shown  in  figure  16,  the  function  approaches  zero 
but  never  crosses  the  axis.  This  sort  of  false  root  will  cause  the  root-finder  to  fail  to 
converge  in  the  specified  number  of  maximum  iterations.  This  problem  is  avoided 
through  the  use  of  accurate  initial  values.  However,  the  false  root  may  approach 
the  legitimate  root  over  the  course  of  the  calculations.  In  some  cases,  the  false  root 
may  even  cross  over  the  legitimate  root.  This  is  remedied  by  ensuring  that  the 
root  finder  will  fail  to  converge  in  the  maximum  number  of  iterations  only  for  the 
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case  of  false  roots.  If  the  maximum  number  of  iterations  is  reached,  a  false  root 
is  assumed  and  a  new  starting  value  will  be  chosen  based  on  the  trend  of  the  last 
two  roots  found  for  that  particular  band.  This  jump  bypasses  the  false  root  and 
facilitate  convergence  on  the  legitimate  root. 
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Figure  15  Plot  of  the  determinant  of  the  30A  Si/SiGe  quantum  well  boundary- 
condition  matrix  vs.  energy. 
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Figure  16  Plot  of  the  determinant  of  the  30A  Si/SiGe  quantum  well  boundary 
condition  matrix  vs.  energy. 

3.3  Code  Generation 

The  computational  tools  used  in  this  research  were  developed  in  FORTRAN 
77.  Several  linear  algebra  tools  were  required  to  perform  tasks  such  as  determin¬ 
ing  the  eigenvalues  and  eigenvectors  of  complex  non-Hermitian  matrices.  Several 
options  available  were  Eigenvalue  Subroutine  Package  (EISPACK),  Linear  alge¬ 
bra  Package  (LINPACK),  and  Linear  Algebra  Package  (LAPACK).  LAPACK  was 
chosen  as  it  was  the  most  modern  and  the  most  recently  updated  of  the  available 
subroutines. 

A  software  package  was  developed  to  compute  bandstructures  for  bulk  ma¬ 
terial  and  semiconductor  quantum  wells.  The  package  consists  of  six  programs 
titled  EFA110,  EFA001,  DISP110,  DISP001,  BANDS001,  and  BANDS110. 
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BANDS001,  and  BANDS110  compute  strained  bulk  material  band  structures  for 
the  [110]  and  [001]  growth  directions  respectively.  EFA110  and  EFA001  compute 
the  determinant  of  the  quantum  well  boundary  condition  matrix  as  a  function  of 
energy  for  any  single  value  of  k||.  DISP110  and  DI SP001  utilize  seed  values  from 
EFA110  and  EFA001  to  determine  the  bandstructure  for  any  azimuthal  direction 
in  the  plane  of  the  quantum  well.  All  programs  and  subroutines  developed  are 
listed  in  Appendix  A  with  a  brief  description  of  each. 
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IV.  Results  and  Conclusions 

This  chapter  presents  the  bulk  material  band  structure  and  the  quantum  well 
bandstructure  for  several  materials.  For  the  quantum  wells,  both  the  [110]  and  the 
[001]  growth  directions  are  evaluated. 

4.1  Bulk  Material  Band  Structures 

The  band  structure  for  alloys  consisting  of  strained  Si/Si7oGe3o  ,  GaAs,  and 
In3oGa7oAs  are  plotted  from  data  files  generated  using  the  code  described  in  the 
previous  chapter.  Bands  for  several  directions  in  k- space  are  evaluated.  The  zone 
boundary  for  these  materials  is  approximately  0.5  Angstrom _1  as  given  by  |  where 
a  is  the  material’s  lattice  constant. 

K  r  L  K  r  x 
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Figure  17  Bulk  material  dispersion  curves  for  [001]  strained  GaAs. 
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Figure  18  Bulk  material  dispersion  curves  for  [001]  strained  In30Ga7oAs. 
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Figure  19  Bulk  material  dispersion  curves  for  [001]  strained  Si7oGe30. 


4-2  Quantum  Well  Band  Structures 

The  band  structure  for  several  types  of  quantum  wells  are  shown  for  Si/Si7oGe30 
,  Al3oGa7oAs/GaAs,  and  GaAs/In3oGa7oAs.  The  following  plots  show  the  evolu¬ 
tion  of  the  bandstructure  under  an  azimuthal  rotation  in  k  — space.  By  noticing  the 
evolution  of  the  band  structure  as  the  azimuthal  angle  progresses,  the  magnitude 
of  the  anisotropy  can  be  observed.  The  directions  are  defined  with  respect  to  the 
2  axis  which  corresponds  to  the  growth  direction,  ky  has  the  components 

kx  =  |ky|cos0 
ky  =  |ky|sin0 
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The  band  structures  in  the  following  plots  are  shown  as  a  function  of  the  azimuthal 
angle  9.  The  final  plots  show  surfaces  representing  the  bandstructure  in  the  positive 
kx  and  ky  quadrant. 
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Figure  20  Dispersion  curves  for  [001]  Al3oGa70As/GaAs  quantum  well  with  L=50 
angstrom.  Theta  ranges  from  0  (dark  line)  to  45  degrees  in  10  degree  increments 
(dashed  lines). 
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Figure  21  Dispersion  curves  for  [110]  Al30Ga7oAs/GaAs  quantum  well  with  L=50 
angstrom.  Theta  ranges  from  0  (dark  lines)  to  90  degrees  in  10  degree  increments 
(dashed  lines). 
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Figure  22  Dispersion  curves  for  [001]  In30Ga7oAs/GaAs  quantum  well  with  L=50 
angstrom.  Theta  ranges  from  0  (dark  line)  to  45  degrees  in  10  degree  increments 
(dashed  lines). 
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Figure  23  Dispersion  curves  for  [110]  In30Ga7oAs/GaAs  quantum  well  with  L=50 
angstrom.  Theta  ranges  from  0  (dark  line)  to  90  degrees  in  10  degree  increments 
(dashed  lines). 
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Figure  24  Dispersion  curves  for  [001]  Si/Si7oGe3o  quantum  well  with  L=30 
angstrom.  Theta  ranges  from  0  (dark  line)  to  45  degrees  in  10  degree  increments 
(dashed  lines). 
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Figure  25  Dispersion  curves  for  [110]  Si/Si7oGe3o  quantum  well  with  L=30 
angstrom.  Theta  ranges  from  0  (dark  line)  to  90  degrees  in  10  degree  increments 
(dashed  lines). 
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Figure  26  Surface  plot  of  the  bandstructure  for  [001]  AboGaroAs/GaAs  quantum 
well  with  L=50  angstrom. 
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4-3  Anisotropy 

Non-parabolic  bandstructure  is  a  result  of  the  periodic  potential  created  by 
the  atomic  arrangement  in  the  crystal  lattice.  The  atomic  arrangement  is  not 
isotropic  and  therefore  the  periodicity  of  the  potential  is  not  isotropic.  In  turn,  the 
bandstructure  is  not  isotropic. 

From  the  plots  show  in  the  previous  section,  a  90  degree  azimuthal  rotation 
in  the  plane  of  the  quantum  well  has  different  results  for  the  [110]  and  the  [001] 
materials.  The  [001]  crystals  have  identical  bandstructures  for  9  =  0  and  9  =  90. 
The  maximum  variance  in  the  bandstructure  from  that  at  either  of  these  values 
of  9  occurs  at  45  degrees.  For  the  [110]  materials,  the  anisotropy  is  even  more 
significant  as  the  bandstructure  continuously  changes  through  the  entire  90  degree 
rotation. 

The  band  structure  in  each  growth  direction  in  each  type  of  quantum  well  for 
each  type  of  material  clearly  shows  an  anisotropy.  It  can  also  be  reasonable  inferred 
that  the  anisotropy  in  the  bandstructure  is  indicative  of  the  behavior  of  the  wave 
functions.  Therefore  calculations  that  require  knowledge  of  the  band  structure 
and/or  the  wavefunctions,  such  as  that  of  absorption  coefficients,  must  consider 
this  anisotropy.  However,  also  observed  from  the  plot  is  that  the  anisotropy  is 
not  significant  for  small  values  of  ky.  This  implies  that  transitions  dictated  by  the 
Fermi  energy  to  be  near  ky  «  0,  can  safely  neglect  the  effects  of  bandstructure 
anisotropy. 

4-4  Summary  and  Recommendations 

A  software  package  has  been  developed  to  evaluate  bandstructures  in  semi¬ 
conductor  quantum  wells  for  materials  grown  in  both  the  [110]  and  the  [001]  di¬ 
rections.  These  tools  have  allowed  analysis  of  the  anisotropy  of  the  quantum  well 
bandstructure  in  the  plane  of  the  well. 
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To  fully  understand  the  magnitude  of  the  anisotropic  effects,  absorption  co¬ 
efficients  themselves  should  be  calculated.  The  code  generated  during  this  research 
provides  everything  necessary  to  perform  the  calculation  of  the  wave  functions  in 
the  well.  From  this,  the  momentum  matrix  elements  and  finally  the  absorption 
coefficients  can  be  calculated.  The  tools  developed  in  this  research  can  also  be 
used  for  the  generation  of  plots  depicting  Fermi  surfaces.  One  final  opportunity 
for  further  research  is  to  include  the  effect  of  valence  band  and  conduction  band 
coupling  by  expanding  the  total  Hamiltonian. 
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Appendix  A.  Subroutine  Description 

The  following  describes  the  programs  and  subroutines  generated  to  perform  the 
necessary  calculations  in  this  research. 

•  EFAxxx  :  This  program  generates  a  data  list  for  plotting  the  determinant 
of  the  quantum  well  boundary  condition  matrix  verses  energy.  The  data  list 
provides  three  records  per  field:  energy,  real  part  of  the  determinant,  and  the 
imaginary  part  of  the  determinant. 

•  OFFSETxxx  :  This  subroutine  determines  the  value  of  the  constants  Vb 
and  Vw. 

•  FACTxxx  :  This  subroutine  utilizes  the  material  constants  to  develop  the 
matrix  elements  for  the  matrices  making  up  the  factored  Hamiltonian,  HO, 
HI,  and  H2. 

•  KZEMATRIXrnr  :  This  subroutine  utilizes  the  factored  Hamiltonian  matri¬ 
ces  to  develop  the  matrix  equation  describing  kz  as  a  function  of  energy. 

•  ORDER  :  This  subroutine  takes  the  eigenvalues  and  eigenvectors  of  the 
kz(E )  matrix  and  selects  the  linearly  independent  solutions.  These  eigenval¬ 
ues  are  then  ordered  as  HH,  LH,  LH,  HH,  SO,  and  SO.  The  corresponding 
eigenvectors  are  ordered  in  columns  to  create  a  matrix  for  use  in  the  boundary 
condition  matrix  equation. 

•  BCMATRIX  :  This  subroutine  utilizes  elements  form  the  factored  Hamil¬ 
tonian  and  eigenvalue  and  eigenvectors  of  the  kz(E )  matrix  to  develop  the 
matrix  equation  describing  the  boundary  conditions  for  the  quantum  well. 

•  DISPxxx  :  This  program  calculates  quantum  well  band  structures.  It  incor¬ 
porates  a  modified  secant  method  root  finder  to  look  for  roots  at  each  step  in 
k  parallel.  The  output  from  this  program  is  a  data  file  containing  the  value 
of  k  parallel  and  the  associated  roots. 
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•  optionsxxx x  :  This  is  the  input  data  file  that  is  read  by  both  EFA  and  DISP. 
The  options  include  material  type,  well  width,  number  of  bands,  and  seed 
values  for  the  root  finder.  Only  material  type  and  well  width  are  used  in 
EFA. 

•  material  :  This  is  a  data  file  specified  in  the  options  file.  It  contains  stoi¬ 
chiometry  values,  Luttinger  coefficients,  deformation  potentials,  average  va¬ 
lence  band  energy,  etc. 
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Appendix  B.  Values  for  Physical  Constants 

The  constants  used  in  the  calculations  were  obtained  from  tabulated  data  in  People 
and  Jackson  (8)  and  Loehr  (7). 


constant 

units 

Si 

Ge 

GaAs 

AlAs 

InAs 

Luttinger  coefficients 

7i 

4.22 

13.4 

7.65 

3.45 

19.67 

72 

0.39 

4.25 

... 

2.41 

0.68 

8.37 

7s 

1.44 

5.69 

3.28 

1.29 

9.29 

lattice  constant 

a 

A 

5.431 

5.658 

5.654 

5.6 

6.0584 

spin-orbit  interaction  term 

A 

eV 

0.044 

0.29 

0.34 

0.28 

0.38 

valence  band  hydrostatic  potential 

<&v 

eV 

2.1 

2.00 

1.16 

2.47 

1.00 

valence  band  uniaxial  deformation  potentials 

Du 

eV 

3.41 

3.32 

2.55 

2.25 

2.70 

D  „/ 

eV 

4.32 

3.81 

3.94 

2.95 

3.12 

elastic  constants 

Cll 

1012  Dyne/cm2 

1.675 

1.315 

1.1879 

1.25 

0.8329 

C12 

1012  Dyne/cm2 

0.65 

0.494 

0.5376 

0.534 

0.4526 

C44 

1012  Dyne/cm2 

m 

0.594 

0.542 

0.396 

average  valence  band  energy 

Ev 

av 

eV 

-7.03 

-6.35 

-6.92 

-7.49 

-6.67 
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Values  for  standard  physical  constants  are  given  below. 


Planck’s  constant 

h 

J  s 

1.05457  x  10"34 

ev  s 

6.58212  x  10~16 

Electron  Mass 

m0 

kg 

9.10939  x  10"31 

eV  (s/A)2 

5.68564  x  10"32 
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